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I. INTRODUCTION 

First search for <f> —> ir°ir°j decay was carried out 
with the ND detector at VEPP-2M e+e" collider in 
1987 J|,§]. In this early experiment the upper limit 
B((f> — > 7r°7r°7) < 10~ 3 was placed. As it was shown later 
by Achasov [p| , study of this decay can provide a unique 
information on the structure of the light scalar /o(980) 
meson. Subsequent studies proved this idea. In 

these works different models of the /o(980)-meson struc- 
ture were considered. The most popular were 2-quark 
model 0, 4-quark MIT-bag model 0, and KK molec- 
ular model Q. 

In 1995 new Spherical Nonmagnetic Detector (SND), 
having better hermeticity, granularity, energy and spatial 
resolution than previous ND detector, started operation 
at VEPP-2M. First indications of the process 



e + e 



^ 7r° 7 . 



(3) 



(1) 



were seen by SND in 1997 (l2|. This preliminary re- 
sult was based on analysis of half of the </>-meson data 
recorded in 1996 - 1997 [jl3|. Present work is based on 
the analysis of full data sample. The — > 7r°7r°7 decay 
was observed for the first time and its branching ratio 
was measured. 



II. DETECTOR AND EXPERIMENT 

The SND detector [Q is designed for the study of 
e + e~-annihilation at center of mass energy about 1 GeV . 
Its main part is a three layer electromagnetic calorimeter, 
consisted of 1630 Nal(Tl) crystals The energy reso- 
lution of the calorimeter for photons can be described 
as <j e /E = 4.2%/ yE{GeV) angular resolu- 

tion is close to 1.5 degrees. The 1996 experiment con- 
sisted of 7 independent runs jl3). In each run the data 
were recorded at 14 different beam energies in the region 
2E Q = (980 - 1050) M eV covering the peak and close 
vicinity of the ^-resonance. The energy spread of the 
beam energy was equal to 0.2 MeV. The total integrated 
luminosity of 3.9 pb^ 1 was collected, corresponding to 
about 8 • 10 6 4> mesons produced. The luminosity de- 
termination was based on events of Bhabha scattering 
and two-photon annihilation detected in the calorimeter. 
Accuracy of normalization is better than 3 %. 



III. DATA ANALYSIS 

Main resonant background to the decay ([!]) comes from 
the decay 



e —> <j) VI ~ > 37r°7 



(2) 



due to the merging of photons and/or loss of photons 
through the openings in the calorimeter. The main 
source of non-resonant background is a process 



The background from the </> — * pir° — * 7r°7r°7 decay is 
small ppq] , nevertheless its amplitude was included into 
simulation of the process (|^). The background from the 
QED 57-annihilation process was estimated and found to 
be negligible. 

In order to suppress the contribution of background 
the events were selected with 5 photons emitted at an- 
gles of more than 27 degrees with respect to the beam 
and without charged particles. Standard SND cuts on 
energy-momentum balance in an event and on photon 
identification function were used: 

• 0.8 < E dep /2E < 1.1 , where E is the beam 
energy, Ed ep - the total energy deposition in the 
calorimeter; 

• p/2Eq < 0.15, where p - is a total momentum of 
all photons in the event; 

• X 2 < 20 - event kinematics is consistent with 7r°7r°7 
hypothesis; 

• C<0; 

The £ parameter was constructed on the basis of the 
logarithmic likelihood function, describing the proba- 
bility for observed transverse profile of electromagnetic 
shower in the calorimeter to be generated by a single 
photon |l^JT^]. This parameter facilitates efficient sep- 
aration of events with well isolated photons from those 
with merged photons or with showers, produced by Kl 
mesons. 

The x 2 is a kinematic fit parameter, describing the 
degree of energy-momentum conservation in the event 
with additional requirement that two 7r° mesons are 
present. During this fit all possible combinations of pho- 
tons (Fig.|l|) in the event were checked in a search for 
invariant masses mi and mi, satisfying the condition: 
y/(m\ - m w ) 2 + (m 2 - m^) 2 < 25 MeV/c 2 . 

The <fi — * KsKl — ► tt°tt°Ki j decay can contribute due 
to nuclear interactions of mesons in the material of 
the calorimeter but after described cuts it is not seen at 
present level of statistics. 

In the energy region of this experiment the invariant 
mass of the pion pair in the process (||) is less than 
700 MeV. In the events satisfying this condition clear 
w(782) peak is seen in distribution (Fig.|J), prov- 
ing the domination of the process (3|) in this kinematic 
region. The m X7 parameter was defined as an invari- 
ant mass of the recoil photon and one of 7r° mesons, 
closest to the w-meson mass. 499 found events with 
750 MeV < m-n-y < 815 MeV were assigned to ujtt" 
class, while 189 events with m I7 outside this interval 
and ttitxt! > 700 MeV were assigned to ir°ir°j class. Sub- 
tracting calculated contribution of the process (||) and 
using estimated probabilities of events misidentification 
for the processes ([!]) and (||) the number of the events of 
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the process @ in the u>ir° class was estimated to be equal 
to 449. The corresponding number of events of the decay 
(|l|) in the 7r°7r°7 class is 164. The background from the 
process (||) was estimated using events of the ujtt class. 
No additional knowledge of the actual production cross 
section of this process was necessary. 

Then, for the events of the ir°ir°j and ojir classes the 
comparison of experimental and simulated distributions 
in ip and 9 angles was done. Here ip is an angle of the 
recoil photon with respect to pion direction in the 7r°7r° 
center of mass reference frame, 6 is an angle between 
recoil photon and the beam. The distribution in 9 for 
7r 7r 7 events with pion pair in a scalar state must be 
proportional to 1 + cos 2 9 and uniform in cosip. The com- 
parison (Fig.||a,c) shows, that in 7r°7r°7 class of events 
pions are actually produced in scalar state. On the con- 
trary, the experimental events of the ujtt class (FigJ^b,d) 
well match the hypothesis of the intermediate ujtt state 
with quite different ip distribution. 

The 7r°7r° invariant mass distribution for the events 
with m^y outside the 750 MeV < m^ 7 < 815 MeV inter- 
val (Fig.Ek) shows significant excess over background at 
large m wn . At m W7r < 600 MeV the sum of background 
contributions dominates. Detection efficiency (Fig.^b) 
for the process (]l|) was determined using simulation of 
the process <p — > Sj — > 7r°7r°7, where S is a scalar state 
with a mass ranging from 300 to 1000 MeV and zero 
width. In addition, this simulation provided information 
on 7r°7r° invariant mass resolution and event misidentifi- 
cation probability as a function of Tn W7t . 

After background subtraction, correction for detection 
efficiency, and correction for misidentification the mass 
spectrum was obtained (Table 1). For masses in the 
(600 — 850) MeV interval the invariant mass resolution 
is equal to 12 MeV , so the 20 MeV bin size was cho- 
sen. At higher masses the resolution improves, reaching 
7.5 MeV at 950 MeV, thus the bin size of 10 MeV was 
used for higher masses. 

To check accuracy of the detection efficiency, obtained 
by simulation, events of the process (0) with 7 photons 
in the final state and reconstructed 3 7r mesons were an- 
alyzed (20) . All other selection criteria were the same as 
in the 7T^r°7 analysis. The number of observed events 
of the process (||) together with the PDG Table value 
for B{<j> rn) = (1-26 ± 0.06) % (UJ and known to- 
tal integrated luminosity provided independent efficiency 
estimation, based on experimental data. Obtained 10% 
difference between experimental and simulated detection 
efficiencies was used as a correction to the 7r°7r°7 detec- 
tion efficiency, obtained by simulation. Such an approach 
minimizes systematic errors corresponding to inaccurate 
simulation of tails of distributions over parameters, used 
for event selection. 



IV. RESULTS 

Summing data from the Table 1, one can estimate 
branching ratio of the decay ([j]) for the mass range 
TTiff-x > 700 MeV 

B{<P -> 7r°7r°7) = (1.00 ± 0.07 ± 0.12) • 10~ 4 , (4) 

and for > 900 MeV 

B{<P -> 7r°7r°7) = (0.50 ± 0.06 ± 0.06) • 10" 4 . (5) 

Here the first error is statistical and the second is sys- 
tematic, which was estimated to be close to 12 %. The 
systematic error is mainly determined by the following 
contributions: 

• the background subtraction error, which decreases 
almost linearly with the increase of invariant mass 
m I7r and is 5 % on average; 

• error in the detection efficiency estimation, which 
increases with and is equal to 8 % on average; 

• systematic error in B((p — > 777) is equal to 5 %. 

The TOtttt invariant mass spectrum (Fig.^J) was fitted 
according to Refs. || and further used for simulation of 
the decay ([!]). As a result the detection efficiency of the 
process ([j]) was estimated (14.72 %) for invariant masses 
within the (600 - 1000) MeV interval (Fig{|a). This ef- 
ficiency was used in the fitting of the c/)-resonance excita- 
tion curve. The visible cross section in each energy point 
&vis{s) was described as a sum of the processes (jlj), (^), 
and (|3|) with efficiency and radiative corrections taken 
into account for each process. The background cross sec- 
tion due to the process (§) was estimated by fitting the 
visible cross section of the events of the ujtt type with 
linear function. The background from the process (||) 
was obtained from simulation. The only free parameter 
of the fit was the <p — > 7r°7r°7 branching ratio, all other 
0— meson parameters were taken from the PDG data [5TJ . 
As a result (Fig.||), the following value was obtained: 



7r°7r°7) = (1.14 ± 0.10 ± 0.12) • 10" 



(6) 



which, in contrast with (Q) and (g), is valid for the whole 
mass spectrum. In the systematic error estimation the 
following considerations were taken into account. In com- 
parison with the results (||) and (|||) the accuracy of nor- 
malization 3 % and efficiency estimation 5 % are higher 
here, the background subtraction error 5 % is the same, 
but an additional systematic error 6 % exists, due to un- 
certainty in extrapolation of the invariant mass spectrum 
into the region < 600 MeV. Smaller systematic un- 
certainty has a ratio of branching ratios: 



7r°7r°7) 



B(4> -> 777) 



(0.90 ±0.08 ±0.07) • 10" 



(7) 
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Assuming that the process (Q) is fully determined by 
/07 mechanism, using the relation B(fg — > 7r + 7r~) = 
2S(/ -> 7r°7r°), and neglecting the decay (f> -> if ^7 [|, 
we can obtain from (|^) 

-8(0 -> /o (980)7) = (3-42 ± 0.30 ± 0.36) • 10~ 4 . (8) 

To monitor stability of experimental results, the number 
of 7r°7r°7 events was checked separately in all 7 experi- 
mental runs. The results for all runs agree well within 
statistical uncertainty. Another important proof of va- 
lidity of the results is the study of the process (0) M , 
which properties are very close to the decay under study. 
The methods of analysis are also very close, and the mea- 
sured value of the decay branching ratio (^) agrees with 
the PDG table value 

V. DISCUSSION 

We would like to emphasize that all results presented 
in the previous section are not based on model assump- 
tions about /o structure. The enhancement at large m vv 
(Fig.[|) is compatible only with large /07 contribution. 
Further analysis is carried out under assumption that the 
decay mechanism is a pure /07 and the mass spectrum 
can be described by the following expression p|,^8[: 

dBrjcp -> ^ 7r° 7 ) = 2m 2 r(/ -> 7r°7r°)r(0 -> / o7 ) 
dm^ n\Df (m T7r )| 2 

where Df{m vv ) is an fo propagator, partial widths 
r(/o — > 7r°7r°) and T((j) — > /07) are functions of m^. In 
the "narrow resonance" approximation [^2| the results of 
the fit of the spectrum are the following: 

m/ = (984 ± 12) MeV, T f = (74 ± 12) MeV 

in agreement with PDG data The coupling constant 
of fo with pions is found to be g , ? 7r+7r _ /A-k — (0.20 ± 
0.03) GeV 2 . 

On the other hand this simple approximation seems 
inadequate for fo meson. In this case width corrections 
might be very large J2^] . Instead of calculating such cor- 
rections we used complete formulas from Refs. ||[l8| for 
the "broad resonance" fit. The formulas take into ac- 
count final width of fo and strong coupling with KK. 
The results of this fit are the following: 

to/ = (971 ±6) MeV, T f {m s ) = (188±||) MeV, 

= 2.10™ GeV 2 , = 0.51±g:g GeV 2 , 



in Figj^ practically coincide, although the "broad reso- 
nance" parameterization looks more realistic. In the lat- 
ter case the value of the coupling constant g 2 KK /^ ob- 
tained from the fit agrees with the predictions of 4-quark 
MIT-bag model (2.3 GeV 2 |Q) as well as the value of 
the branching ratio (j^). The corresponding predictions 
of the KK molecular model [p|]T8]| and the 2-quark model 
H are about 5 times lower. 

Of course, besides fo other intermediate states may 
contribute to the decay under study, and most proba- 
bly the heavy and broad c-state |l8| ]. This contribution 
can decrease values of fo coupling constants, but in any 
case such a state must produce nearly flat 7r°7r° invariant 
mass spectrum, not masking /o-meson resonance signal 
in the mass region > 900 MeV. Theoretical esti- 

mations of branching ratios are also imprecise depend- 
ing on a set of not well defined coupling constants and 
may vary within a factor of two Jl8[ . Thus, measured 
<j> — » /q7 branching ratio, which is 5 times higher than 
2-quark and KK model predictions can be considered as 
another indication of significant 4-quark MIT-bag part in 
the fo meson not excluding participation of other quark 
configurations. To determine contributions of all possible 
transition mechanisms and all quark configurations in fo 
one has to perform simultaneous analysis of all fo data. 

VI. CONCLUSION 

In the SND experiment the <f> — > 7r°7r°7 radiative decay 
was observed for the first time and its branching ratio 
was measured. Observed enhancement at high to^ in 
the invariant mass spectrum together with angular dis- 
tributions agreeing with scalar intermediate 7r°7r° state 
show, that /o (980)7 transition mechanism dominates in 
this decay. 

Measured fo parameters are strongly model dependent 
due to the fact that fo meson is broad and its mass is 
close to KK threshold, making simple Breit-Wigner de- 
scription of the resonance inadequate. 

Observed high decay probability ~ 10~ 4 gives evi- 
dence, that fo meson contains significant 4-quark com- 
ponent but does not exclude participation of other quark 
configurations. 
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9 2 f K+K -/9% + n- =4.1 ±0.9. 

Comparison of the results of the fits shows strong model 
dependence of /o-meson parameters. Both fitting curves 
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FIG. 1. Invariant mass spectrum of photon pairs for se- 
lected 5-gamma events. 
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FIG. 2. Distribution over it 7 invariant mass for 
ttitttv < 700 MeV. Points - data, histogram - simulation, 
shaded histogram - sum of simulated contributions from 
4> ~> VI an d 4> — * 7i" f V°7 decays, arrows - selection of the 
our class. 
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FIG. 3. a, b - cosine of ip, the angle between directions 
of 7T° and recoil 7 in the rest frame of 7r°7r° system; c, d - 
distributions in 9, angle of the recoil 7 with respect to the 
beam. Points - data, histogram - simulation. 
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FIG. 4. a - invariant mass distribution of 7r°7r° pairs for 
selected 7r°7r°7 events without acceptance corrections. His- 
togram - data, shaded histogram - estimated background 
contribution from e + e~ — » con and <j> — > 777; b - detection 
efficiency for 7r°7r°7 events. 




FIG. 6. Energy dependence of the visible e e 
cross section. Points - data, solid line - fit, dotted line - 
estimated background contribution from e + e~ — > uitt and 

<f> — > 777. 
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FIG. 5. The measured n n invariant mass spectrum. 
Background is subtracted and efficiency corrections applied. 
Points - data, solid line - the result of the "broad resonance" 
fit, dashed line - the result of the "narrow resonance" fit. 



TABLE I. The experimental mass spectrum for the 
<f> — » 7r°7r°7 decay after background subtraction and accep- 
tance corrections. Only statistical error are shown. 



10 7 (MeV 



(MeV) 


dBr{4> — *n° -K ~y) 


600-620 


0.61 ± 0.84 


620-640 


1.26 ±0.89 


640-660 


0.59 ± 0.68 


660-680 


1.66 ± 0.84 


680-700 


1.62 ±0.78 


700-720 


0.81 ±0.61 


720-740 


2.86 ± 0.94 


740-760 


1.63 ± 0.76 


760-780 


1.73 ± 0.80 


780-800 


2.69 ± 0.96 


800-820 


2.50 ± 0.94 


820-840 


2.25 ± 0.90 


840-860 


4.58 ± 1.21 


860-880 


2.32 ± 0.91 


880-900 


3.80 ± 1.11 


900-920 


3.59 ± 1.14 


920-940 


3.74 ± 1.17 


940-950 


8.41 ± 2.47 


950-960 


8.93 ± 2.62 


960-970 


6.50 ± 2.34 


970-980 


3.74 ± 1.93 


980-990 


6.20 ±2.57 


990-1000 


1.60 ± 1.60 
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